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ABSTRACT 

While it is incontrovertible that the inner Galaxy contains a bar, its structure near the Galactic 
plane has remained uncertain, where extinction from intervening dust is greatest. We investigate 
here the Galactic bar outside the bulge, the long bar, using red clump giant (RCG) stars from 
UKIDSS, 2MASS, VVV, and GLIMPSE. We match and combine these surveys to investigate 
a wide area in latitude and longitude, \b\ ^9° and |/| ^ 40°. We find: (i) The bar extends to 
I ^ 25° at jZ?! ^ 5° from the Galactic plane, and to I ^ 30° at lower latitudes, (ii) The long bar 
has an angle to the line-of-sight in the range (28 — 33)°, consistent with studies of the bulge 
at |/| < 10°. (iii) The scale-height of RCG stars smoothly transitions from the bulge to the 
thinner long bar. (iv) There is evidence for two scale heights in the long bar. We find a ^ 180pc 
thin bar component reminiscent of the old thin disk near the sun, and a ^ 45 pc super-thin 
bar component which exists predominantly towards the bar end. (v) Constructing parametric 
models for the RC magnitude distributions, we find a bar half length of 5.0 ± 0.2kpc for the 
2-component bar, and 4.6 ± 0.3kpc for the thin bar component alone. We conclude that the 
Milky Way contains a central box/peanut bulge which is the vertical extension of a longer, 
flatter bar, similar as seen in both external galaxies and N-body models. 
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1 INTRODUCTION 


Gas kinematics and near-infrared (NIR) photometry {e.g. |Binney| 
|et al.|199l]|Weiland et al.|1994) have shown that the Galactic bulge 
at longitudes |/| < 10° contains a bar-like structure. These results 
have since been confirmed, and the barred bulge has been character¬ 
ised and mapped with increasing accuracy. Using star counts of red 
clump giants (RCGs), [Stanek et al.| ( |1994| [19971 1 reached a similar 
conclusion that the bulge was barred, a result since confirmed with 
increasing confidence and accuracy ( |Babusiaux & Gilmore|2005] 
[Rattenbury et al.|2007| l. The studies using RCGs agree with results 
from brighter giant star counts ( [Lopez-Corredoira et al.|2005]l and 
the COBE integrated NIR emission (|Dwek et al.|19^ Binney et^ 


|1997|[Freudenreich|1998|[Bissantz & Gerhard|2002| ). Most authors 

have found a bar angle, defined as the angle of the major axis of the 
bar to the sun line-of-sight, in the range 20 — 30°. Recent interest 
was stimulated by the discovery that close to the bulge minor axis 
at \b\ >5° the RCGs have two distinct magnitudes (MeWilliam & 


|Zoccali|2010[|Nataf et al.|2010] t. This is because these lines-of-sight 
pass through both arms of an X-shaped structure which is charac¬ 
teristic of boxy/peanut (B/P) bulges in barred galaxies ( |McWilliam| 
|& Zoccali|2010|[Ness et al.|2012| ). These B/P bulges arise naturally 
in N-body simulations of disk galaxies {e.g. |Athanassoula|2005l 
[Martinez- Valpuesta et al.|2006]l and are comm on in external galax- 
ies ( [Bureau et al.|2006[|Laurikainen et al.|2014| t. This has stimulated 
further investigation of the full three dimensional structure of the 
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bulge beginning with [Saito et al.[ ( [20TT| l using 2MASS data. Vista 
Variables in the Via Lactea (VVV) has provided significantly deeper 
and more complete data and this was recently used to map the B/P 
bulge non-parametrically in three dimensions by [Wegg & Gerhard| 

It has also been suggested in NIR star counts using in-plane 
data at 10° < / < 30° that there is a less vertically extended 4 kpc 
bar at ^45° to the line of sight ( [Hammersley et al.|1994[ l. This 
component has been termed the long bar. The existence of the 
long bar was confirmed using RCGs by [Hammersley et W. { 


poooli 


and subsequently with increasingly more powerful NIR ( Lopez-[ 
[Corredoira et al.[200^[Cabrera-Lavers et al.[2007([2008[ t and longer 


wavelength GLIMPSE data ( [Benjamin et al.[2005[[Zasowski[2012[ ). 

Understanding the nature of the long bar and its structure, 
amplitude, length and pattern speed is of great importance for many 
areas of Milky Way study. It influences for example the disk outside 
the bar ( [Minchev & Famaey[[2010[ l, the kinematics in the solar 
neig hbourhood ([Dehnen[2000| l, and the observed non-circular gas 
flow ( [Bissantz et al.[2003k 

One important unresolved issue is the relationship between 
the bar in the Galactic bulge at |/| < 10°, and the long bar outside 
it. Studies of the long bar at / > 10° have often found a larger bar 
angle (typically 45°, beginning with [Hammersley et al.[1994] l in 
comparison to the three-dimensional bulge (typically ^25 — 30° e.g. 
(27 ±2)°: [Wegg & Gerhard 20T3] ). It has therefore been suggested 
that the Galaxy contains two bars, with the central three-dimensional 
bar not aligned to the long bar. However an in-plane long bar mis¬ 
aligned with a triaxial three-dimensional bar is difficult to reconcile 
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Figure 1. In the top figure we show the surveys used in this study. We use, in order of preference, VVV in red, UKIDSS in green, and 2MASS in blue. Grey 
regions are those without data of sufficient depth i.e. close to the plane without VVV or UKIDSS data where 2MASS is insufficient. In the lower figure we show 
the surface density of stars in the K^-band in the extinction-corrected magnitude range 12.25 < Kq < 12.75. Asymmetric number counts in / close to the plane 
are a result of non-axisymmetry due to the long bar. The star counts are smoothed with a Gaussian kernel of cr = 0.1°. Extinction is corrected using theH — K 
colour excess as in |equation Kq = -\-Mk^rc) and data outside the colour bar range are plotted at its limit. 


dynamically; the suggested length ratio is low and therefore the 
mutual torques are strong. It has instead been suggested that rather 
than two distinct bars, the long bar is the in-plane extension of the 
central three-dimensional boxy/peanut structure structure ( [Martinez^ 
[Valpuesta & Gerhard|2011|[Romero-G6mez et al.|2011| t. One of the 
motivations for this study is to help resolve this controversy. 

Throughout we use the terminology that the bar outside the 
bulge region at |/| > 10° is the long bar, regardless of the details of 
thickness, bar angle, or alignment with the barred bulge. 

Our primary indicator of bar structure are RCGs which are core 
helium burning stars and provide an approximate standard candle 
( [Stanek et al.|1994| t. We combine several surveys to have the widest 
view and the greatest possible scope on the bar density distribution. 
In the .K^-band we use data from (i) the United Kingdom Infrared 
Deep Sky Survey (UKIDSS) Galactic Plane Survey (GPS, [Luc^ 
|et al. poo’s] !, (ii) the VVV survey ( |Saito et al.|2012| l and, (iii) to 
extend the study further from the galactic plane than previous studies, 
we augment this with 2MASS data ( |Skrutskie et al.|[2006] l. We 
homogenise the analysis of the surveys using a common photometric 
system and identify RCGs statistically in magnitude distributions 
rather than in colour-magnitude diagrams since this has worked well 
in the bulge (e.g. |Nataf et al.|2013[[Wegg & GerhardpOlSj l. 

We verify our results where possible using data at 3.6qm and 
4.5qm, which is significantly less affected by extinction, taken from 
the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) 
survey on the Spitzer space telescope ( [Benjamin et al.|20^5] l. Be¬ 
cause this data only covers \b\ < 1° we use the K^-band as our 
primary data, but the GLIMPSE data remains very important for 
cross checks, particularly of dust extinction. 

This work is organised as follows: In [section^ we describe 
the data and construction of magnitude distributions for the stars 
in bulge and bar fields. In [section"^ we fit the red clump stars in 
these magnitude distributions and discuss the features of these fits. 


In |section'4| we examine the vertical structure of the fitted red clump 
stars in longitude slices, and in |seetion~5] we derive densities whieh 
fit and best fit the observed magnitude distributions. We discuss our 
results and plaee them in context in |section~6| and finally conclude 
in lsection~71 


2 MAGNITUDE DISTRIBUTION CONSTRUCTION 

A number of steps are required to combine the surveys (UKIDSS, 
VVV, 2MASS, GLIMPSE) and bands (//, 3.6|am, 4.5|am) and 

construct consistent magnitude distributions: The surveys must be 
transformed to the same photometric system, extinction corrected, 
and to compare bands and convert to distances we require the char¬ 
acteristic magnitudes and colours of RCGs. 

2.1 H and K^-band data 

The first step in construction of the magnitude distributions is to 
transform all the surveys to the same photometric system. We choose 
to convert the UKIDSS and VVV surveys to the 2MASS system 
using the methods and transformations described in appendix |A| 

2.1.1 Extinction Correction 

Extinction is then corrected for on a star-by-star basis assuming 
that all stars are red-clump giants (RCGs). We primarily work with 
the RCG K^-band distance modulus, fix, where we calculate the 
K^-band extinction from the H — Kg reddening: 

Extinction Correction 

/■ ’S 

Av 

HK=Ks- [(tf-g,)-(tf-g,)Rc] -M^„RC , (1) 

Reddening 
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where {H — Ks)^c is the intrinsic H — Kg colour of RCGs, is 

Ak 

the absolute ^^-band magnitude of RCGs, and ^ is a constant 

which depends on the extinction law. 

The method of correcting for extinction on a star-by-star basis 
in [equation (l)| is similar to other studies with data close to the 
Galactic plane {e.g. [Babusiaux & Gilmore|2005[|Cabrera-Lavers| 
jet al.|200^ . In the Galactic bulge most of the dust lies in a fore¬ 
ground screen making extinction correction using extinction maps 
possible. In-plane this is not true and necessitates the method of 
extinction correction in [equation (1)| 

In the near-infrared bands used in this work, the vast majority 
of the stars have colours close to that produced by a Rayleigh-Jeans 
spectrum, therefore this star-by-star extinction correction using the 
H — Kg is relatively accurate {e.g. the NICE method, p^ada et ah] 
[1994[ ). Moreover, we are concerned with RCGs in this work, which 
as well as being good standard candles naturally have very similar 
colours. For easy comparison between bands we work with the 
distance modulus, under the assumption that the stars are RCGs. 
Other stars form a background which we do not study. We use 
H — Kg io correct rather than J — Kg because the //-band has lower 
extinction and therefore remains deeper in high extinction regions. 
In addition H — Kg more constant across stellar types than J — Kg 
improving the extinction correction ( [Majewski et al.[20TT| l. 

To utilise [elation (l)[ we require two additional ingredients: a 
value for and the magnitude of RCGs for (// — Kg)Rc and 

Mk,,rc- 

To estimate ^ we use the value = 1.73 meas¬ 

ured by [Nishiyama et al.[ ( [2009[ l using red clump stars close to 
the galactic plane towards the galactic center. This corresponds to 

Ak 

E{H-K ) ~ principle the //-band could be used as an ad¬ 

ditional confirmation of our results, however the adopted A///A/^^ 
results in a //-band equivalent to [equation (l)[ for jin which is com¬ 
pletely degenerate with fiK- Instead we use the GLIMPSE survey to 
confirm results where possible because the longer wavelengths are 
less susceptible to dust extinction. 

We estimate the absolute magnitudes of red clump stars in 
a two step process. First we estimate colours of red clump stars 
using the Padova isochrones. Constructing the luminosity function 
from the solar metallicity lOGyr isochrone with a Kroupa IMF and 
fitting a gaussian to the red clump results in red clump colour of 
H — Kg = 0.09. This colour (and the GLIMPSE colours utilised 
later) are largely independent of age and metallicity, changing by 
less than 0.03 mag when varying age from IGyr to 15Gyr and 
metallicity from [Fe/H] = —0.7 to [Fe/H] =0.17. Finally we use 
the calibration Mk^ rq = —1.72 which results in bulge red clump 
stars being approximately symmetric about Rq ~ 8.3kpc ([Wegg &[ 
[Gerhard[2013] l. 


2.1.2 Catalogue Construction 

To utilise the star-by-star extinction correction methods we require 
band matched catalogs to calculate the reddening of each individual 
star. UKIDSS and 2MASS provide band matched catalogues while 
for the VVV survey we match between the H and /T^-band DR2 
catalogs using a radius of 1". 

In combining the surveys we preferentially use VVV data when 
available, falling back to UKIDSS and finally 2MASS. The resultant 
footprint of the surveys is shown in the upper panel of Fig. In 
the lower panel of Fig. we show the surface density of stars in the 
extinction corrected K^-band of our combined data. The long bar 
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Figure 2. The difference in K^-band extinction calculated using two methods. 
Ak,{EHK) is calculated using the H — Kg reddening from [equation (1)[ 
Ak, (rice) is the K^-band extinction calculated by converting the H —[4.5 ji] 
reddening from [equation (2)[ to A^^^. All stars with 12.5 < jik < 14 are 
included. Contours are equally spaced in density at 10, 20... 90% of the peak 
density. Error bars represent the sigma clipped mean and standard deviation 
binned as a function of A^^^ (RICE). 



Figure 3. Example histograms of distance modulus, jd, calculated using 
equations!^ and [^ assuming that all stars are red clump stars. We show 
K^-band (black), 3.6pm (red) and 4.5pm (blue). We also show the fits to 
the histograms made using [equation (3)[ as the solid lines. For legibility the 
3.6pm, and 4.5pm bands are offset by logiV = 0.5 and 1 with respect to the 
axis. This example field has its center at / = 18.5° b = 0.9° and has size 
A/ = 1°, AZ7 = 0.3°. 

corresponds to the non-axisymetric in-plane enhancement of stars 
at positive longitudes. 

2.2 GLIMPSE data 

The longer wavelength GLIMPSE has lower dust extinction and we 
use it to cross check the K^-band data where possible. We use the 
H — [4.5fl] colour to correct for extinction corresponding the RICE 
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Figure 4. Parameters of the Gaussian fits to red clump stars in the ^j-band using equation In the top panel we show the number of red clump stars N^c- In the 
middle panel we show the distance to the fitted peak of the red clump stars, Jpc (calculated from /ipc)- In the lower panel we show the dispersion in magnitudes of 
the fit to the red clump, o^c. Grey regions are those where no well localised fit to the RCG distribution was possible, either because the surveys where not deep 
enough in that field, or the fitted dispersion was too large to represent RCGs towards the Galactic centre. 


method ( [Majewski et al.|20lT| l: 


M3.6m = [3-6jU] - 




-(/f-[4.5p])Hc]-M3*,RC 
4.44. = [4,5;.]) 

-(//-[4.5At])Rc]-M4.5^,RC 


( 2 ) 


where we adopt 


^3.6/f 


= 0.31 and 


^4.5;U 


= 0.23 cor- 


E{H-[A. 5 ii]) ^ 

responding to the measurements of |Zasowski et aL| p009| l in 
the range 15° < / < 20° together with Ah/Ak, = 1.73. Measur¬ 
ing the colour of RCGs using the Padova isochrones as in |Sec-| 


tion 2.1.1 results in [3.6] —Ks = —0.05, [4.5] —Ks = 0.04 and there¬ 


fore M3.6 ^^rc = -1.77 andM4.5^^Rc = -1.68. 

To utilise the H — [4.5fl] colour for extinction correction we 
must band match the GLIMPSE data to surveys providing //-band 
measurements. We have band matched to UKIDSS and VVV H- 
band data, using a radius of 1", while GLIMPSE already provides 
matches to 2MASS. We utilise the deeper UKIDSS and VVV H- 
band data over 2MASS to correct for extinction in GLIMPSE when 
possible. This is particularly important close to the galactic plane 
where otherwise 2MASS would limit the depth of the GLIMPSE 
data to brighter than the clump. 

In |Fig. 2| we compare the K^-band extinction calculated using 
the RICE method ( |equation~(^ ► to the value determined from the 
H — Ks reddening used by [equation (T)| We find good agreement 
between the measurements using the two methods. The mean differ¬ 
ence is less than 0.05 mag at all extinctions. The standard deviation 
in the difference between the estimates is typically less than 0.1 mag 
at magnitudes corresponding to long bar red clump stars. The ma¬ 
jority of this dispersion is due to the photometric error in the 4.5 qm 


band which contributes 0.07 mag at low extinction. We have per¬ 
formed the same comparison as a function of galactic longitude 
and latitude and found no significant differences. This confirms 
that the extinction law does not vary significantly across the region 
considered, and that the different surveys are consistent. 

Finally we construct magnitude distributions in fields covering 
|/| < 40°, \b\ < 9°. We use two grids: close to the galactic plane 
at IZ?I < 1.2° we use fields of size A/ = 1°, AZ? = 0.3°. Away from 
the galactic plane at \b\ > 1.2°, where the number count surface 
densities are lower, we use A/ = 2°, AZ? = 0.6°. 

In |Fig. 3| we show, for an example field, histograms calculated 
using equations [^and]^ In this field RCG stars in the bar are visible 
as the ‘bump’ of stars at q 13.5 above the smooth background of 
non-RCG bar stars. 


3 GAUSSIAN RED CLUMP FITS 


To the magnitude distribution in each field we attempt to fit a Gaus¬ 
sian representing the red clump together with a exponential back¬ 
ground { |Stanek|1995[|Nataf et al.|201^ : 


V(Ai) = 


-^RC 

OrCV^^ 


exp 


1 / M ~Mrc 

2 V tJRc 


+Aexp5(q-qRc) . 

( 3 ) 


We assume that the error is the Poisson error on the number in each 
magnitude bin and fit this equation using Each fit is performed 
using an Markov Chain Monte Carlo (MCMC) of 20,000 steps and 
flat priors, starting from the maximum likelihood position. 

In regions covered by UKDISS or 2MASS we fit over the 
range 11.5 < jj^k < 15.5, in the VVV bulge region we use the 
range 12 < < 16, and in the VVV disk region we use the range 
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Figure 5. In upper half of the figure we show selected slices of the Milky Way bar viewed from above using the ^^-band. In the lower half of the figure we 
show the two component model constructed in |section 5] observed in the same manner. Data are the raw histograms constructed using |equation (l)| with the fitted 
exponential representing non-red clump stars subtracted. The data is divided by the volume of each bin to approximately convert the histograms to densities in 
kpc“^. In the top row of plots we show these densities on a log scale. The lower rows are the same data, but each line-of-sight is normalised to its peak. On the 
0.75 < 1.05 slice we also plot the red clump distances found by |Nishiyama et al!](2009) in blue, and |Gonzalez et al.]j2011) in black, using our adopted 

= —1.72. To reduce noise the histograms are smoothed using boxcar smoothing with a width of 0.15 mag, corresponding to three histogram bins. Red 
lines show longitudes of I = —30, -20,-10,10,20,30° and lines at a = 27° and 45° to the sun line-of-sight with half length 4.7kpc. 


12 < }1k <16.5. These limits were determined by visual inspection. 
The GLIMPSE jU 3 . 6 jU and jUq 5 ^ data is shallower and in this case 
we extend the fit to 10.5 at the bright end, while the faint end is 
empirically determined in each field: In each field we estimate the 
distance modulus at which the counts drops to 50% of the peak due 
to incompleteness and fit to a maximum of 0.75 mag brighter than 
this. Examples of these fits are shown in Fig.[^ 

The resultant fitted parameters, A^rc, Mrc^ and Orc, for the Kg- 
band are shown in Fig.|^ In this plot we only consider the fit as good 
if the red clump is detected at > 3(7 i.e. N^c is at least three times 
larger than its error, and the background slope 5 > 0.55. The second 
criterion is equivalent to rejecting fields which are significantly 
incomplete at the faint magnitude limit since in this case the slope 
of the fitted background is reduced. We also require that the best 
fitting red clump magnitude lies within the range of magnitudes 
fitted, and that the dispersion of the fitted red clump, cjrc < 0.6 
since visually fits that do not meet this criterion appear spurious. 
Finally we exclude all fields at / < —10° without VVV data since 


2 MASS data is not deep enough the detect the clump in this region 
and the few fits which passed the previous criteria were visually 
rejected. We show in |Fig. 6 | the histogram of fitted values of B for 
fields that pass the other selection criteria. 

Several features of the fits in Fig. are worth noting: (i) the 
red clump is detected and fitted up to |Z?| 5° at / > 10° out to 

I ^ 20°, outside the traditionally defined bulge region. In this region 
the fitted distance at the same / is close to the distance fitted in the 
Galactic plane, with an increased dispersion, (ii) Fits no longer pass 
our criteria for a well fitted red clump at / > 30°. This is similar in 
longitude to that identified as the end of the long bar in the literature 
{e.g. |Lopez-Corredoira et al.|2006| Cabrera-Lavers et al.|2007|l. We 
address this in greater detail in |subsection 5.5| and |subsection 6.3| 
(iii) The feature visible from / = —20° to / = —25° in the surface 
density of stars in |Fig. l| is not visible here. This is either because it 
is composed exclusively of young stars which have not yet evolved 
to become clump stars, or is too extended along the line-of-sight to 
give a well localised red clump, (iv) Along the minor axis of the 
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Figure 6. The fitted slope of the background of stars which are not RCGs 
defined though [equation (3^ The grey region are fields which are sub¬ 
sequently rejected for having B < 0.55 and therefore being significantly 
incomplete. The majority of these are crowded central bulge fields. If all 
stars are giants then because of the power law slope of the luminosity func¬ 
tion on the giant branch we would expect this parameter to lie in the range 
0-6 < 5 < 0.78 jMendez et al.T2002) . 


bulge at \b\ >5° the fitted dispersion is larger than off-axis. This is 
a result of using one Gaussian to fit the wider split red clump ( jNatafj 
jet aL|2010| [Me William & Zoccali|2010| l in this region, (v) In the 
Galactic plane at / > 10° the scatter in fitted distance and dispersion 
are considerably smaller than in jCabrera-Lavers et al.| ( |200^ . 

In Fig. we show the raw histograms but with the exponential 
background of non-red clump stars subtracted Le. we subtract the 
fitted A exp5(jU — jUrc) from each histogram. This background sub¬ 
tracted histogram for each line-of-sight is then plotted as if viewed 
from the north Galactic pole. In addition we divide the number 
in each histogram bin by its volume, where the volume of each 
equally spaced bin in magnitude increases with distance like oc . 
This is equivalent to plotting the density assuming that red clump 
stars are perfect standard candles. From these plots it is already 
evident that the bar more closely matches the bar angle shown in 
the figure of a = 27° (g.g. jWegg & Gerhard|2013| ) than ^ 45° {e.g. 
jCabrera-Lavers et al.|2008^ , although the data at low latitudes shows 
a possible curvature towards the end of the bar in a leading sense 
similar to that suggested by [Martinez- Valpuesta & Gerhard| ( |20TT] ). 

We have performed the same analysis with the 3.6qm and 
4.5 qm bands. The equivalent plots to Figs. andare shown fig¬ 
ures |BT|B^in|Appendix^ The GLIMPSE data for b > 0.45° are 
consistent with the .^f^-band data. In particular, they independently 
support a bar angle near 27°. Ai b < 0.45° GLIMPSE does not 
provide sufficient completeness for RCGs for a good comparison. 

In summary we find that the long bar matches a bar angle of 
27° more closely than 45° in both the .^f^-band and the GLIMPSE 
data, and that the long bar is detected to \b\ =5° outside the bulge 
at/<20°. 


4 LONG BAR VERTICAL PROFILE: LONGITUDE 
SLICE FITS 

We estimate the scale height by performing fits to the number of 
red clump stars as a function of latitude. We treat each longitude 
slice independently and fit the N^c found in from equation ([^ to an 
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Figure 7. The vertical profile of the surface density of red clump stars 
identified in the ^^-band in several longitude slices. The blue curve is the 
best fitting single exponential to the .fiT^-band, while the best fitting double 
exponential is the red curve. Error bars show the statistical error resulting 
from fitting [equation (3)[ to each field. 


exponential: 


Vrc(^) = 


bl ) ’ 


(4) 


where bo is an offset and bi the exponential scale height. In fields 
from 15° < / < 25° we find that a single exponential is a poor fit to 
our data, which extends to higher latitudes than previous investiga¬ 
tions. We therefore additionally fit a double-exponential: 


AT / 7 \ ^RC,A 

NRc(b) = ^r^exp 

2/?i,a 


\b-bo\ \ 

blA ) 


^RC,5 

2bi^B 


exp 


\b-bo\ \ 

blA ) 


(5) 

Examples of these fits for the long bar region, where we find 
evidence for two scale heights, are shown in figures and In 

H we show the .^f^-band data and the fits to these profiles. In 
we show the GLIMPSE data but with the fits to the TTy-band 
overplotted. The GLIMPSE data do not extend sufficiently far from 
the galactic plane to verify the double exponential fits, however the 
agreement in the region close to the galactic plane is an important 
check since this is the region where the extinction is greatest and 
the longer wavelength GLIMPSE data is affected less. 

The resultant scale heights of these fits are shown in [Fig. 9[ The 
black points show the exponential scale height of RCGs. The scale 
height increases from the bulge minor axis to a peak near the end 
of the Box/Peanut region, before continuously decreasing through 
the long bar region to < 1° near the bar end. In the long bar region 
where we find that a double-exponential fits the data better the scale 
heights are typically 0.5° and 2°. At / = 20° the bar lies at a 
distance of ^ 5.2kpc and these therefore correspond to scale heights 
of 45pc and 180pc. The 180pc scale-height component is 
similar in scale height to the thin disk in the solar neighbourhood 
and we therefore analogously refer to it as the thin bar. This is 
appropriate since the scale height of thin disks in edge-on external 
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Figure 8. The vertical profile of the surface density of red clump stars in 
several longitude slices in GLIMPSE data compared to the profiles fitted 
to the /^ 5 -band. In black we show the surface density of red clump stars 
identified in the 3.6p.m band and in green those identified in the 4.5p.m band. 
The blue curve is the best fitting single exponential to the ^^-band, while the 
best fitting double exponential is the red curve. The curves have not been 
renormalised. Error bars show the statistical errors resulting from fitting 
jeguation (^ to each field. 



Figure 9. Scale height in degrees in each longitude slice measured by fitting 
the single and double-exponential models to N^c from equation 1^. In black 
we show the results of fitting one exponential, while in red and orange we 
show the results of fitting a double exponential. 


galaxies varies weakly with position fvan der Kruit, P. C. & Searl^ 
|1981[|de Grijs & Peletier|1997) . We refer to the significantly thinner 
^ 50pc component as the super-thin bar, both due to its small scale 
height and in analogy to super-thin components found in edge-on 
external galaxies ( [Schechtman-Rook & Bershady|2013| l. 

We find that the fitted profiles are not symmetric about b = Q 
but instead require an offset. In [Fig. 10| we show the fitted offset 



Figure 10. Longitude offset in degrees in each longitude slice measured by 
fitting the single and double-exponential models to N^c from equation a In 
black we show the results of fitting one exponential, while in red we show 
the results of fitting a double exponential. In blue we show the expected 
offset for a geometrically thin bar with bar angle 27°, the Sun lying 25 pc 
above the Galactic plane, and Sgr A*in the physical Galactic mid-plane at 
b = -0.046° ^Goodman et al.|2014) . 


which is typically ~ —0.1°, corresponding to an offset of 14pc 
at 6kpc. For comparison, we can predict the latitude offset of the 
physical Galactic mid-plane as a function of distance and longitude, 
with two assumptions: that the Sun lies 25pc above the galactic 
plane ( |Mafz- Apellaniz|2QQl[[^en et al.|2001[|Juric et al.|2008| l 
and that Sgr Allies in the Galactic plane with b = —0.046°, so that 
the physical Galactic mid-plane is slightly tilted with respect to the 
b = 0 plane (Fig. 2 of [Goodman et al.|2014| l. We show in blue the 
resultant offset if the bar is geometrically thin along the line-of-sight 
with bar angle 27°. This simple prediction agrees well for / < 15°. 
It slightly overestimates the offset at / > 15° although typically by 
< 0.05°. This could be explained by the apparent bar density peak 
lying behind the geometrically thin prediction for these longitudes, 
as is also seen in jFig. 5| We conclude that the long bar, and barred 
bulge is consistent with lying in the physical mid-plane of the galaxy 
to within approximately 0.05° in latitude, or around 5 pc for these 
distances. This could equally be seen as evidence supporting the 
assumptions that Sgr Allies in the Galactic mid-plane, and the sun 
lies 25 pc above the mid-plane. 

In Fig.[TT]we show the number of stars integrated over Z? as a 
function of longitude and the integrated number of stars converted 
to a number of stars along the bar major axis, assuming a thin 
bar with angle a = 27°. While the errors are large, the number of 
stars in the thin bar component in these plots decreases along the 
bar in a manor consistent with an exponential. Fitting the lower 
panel of [Fig, ll] where we assume that the bar is thin and lies at 
27° to the sun line-of-sight gives the exponential scale length of 
this component to be 1.5kpc, somewhat shorter than the thin disk 
scale length of 2 — 3kpc {e.g. 2.5kpc: [Binney et al.|1997| 2.4kpc: 
Bissantz & Gerhard|2002| 2.6kpc: [Jurid et al.|2008| 2.15^c: |Bov3^ 
& Rix||2Q13j l. The number of stars in the super-thin component 
however increases outwards. We therefore propose that the short 
scale height component corresponds to stars formed more recently 
predominantly towards the bar end, that have not experienced a 
large amount of vertical heating since birth. For reference, to have 
formed RCGs the stars must still be > 500 Myr old but galaxies with 
ongoing star formation have a strong bias to younger RCGs with 
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Figure 11. In the top plot we show the integrated number of RCGs in the 
bar as a function of longitude, Zrc, defined through equatio ns and 
In black we show the result from fitting a single exponential jequation (4^ . 
In red and orange we show the results from fitting a a double exponential 
^equation (5^ in the region where this is a significantly better fit. In the lower 
panel we show the same results with longitude converted to distance along 
the major axis of the bar assuming that it is geometrically thin along the 
line-of-sight and lies at bar angle a = 27°. 


a peak of age ^ 1 Gyr ( [Girardi & Salaris|2001[ [Salaris & Girardi] 
|2QQ2| |. The larger scale height component would then correspond 
to old red clump stars analogous to old thin disk stars in the solar 
neighbourhood. We discuss this scenario further in |section~6| after 
our more detailed modelling. 

We also note that both the scale height and the total num¬ 
ber of red clump stars transitions smoothly from the bulge region 
(|/| < 10°) to the long bar region (10 < / < 30°), providing further 
evidence that the two structures are not distinct. 


5 BEST FITTING PARAMETRIC DENSITY MODELS 

In this section we construct bar/bulge models that match the mag¬ 
nitude distributions and describe their key features such as the result¬ 
ant bar angle, bar length and bar mass. Our approach is to model the 
stellar density, and convolve this with suitable luminosity functions 
to produce model magnitude distributions. The stellar density is 
then adjusted until the model magnitude distributions match the 
extinction corrected magnitude distributions. We again focus on the 


^^-band because of the wider and deeper coverage than GLIMPSE. 
We fit the magnitude distributions over the range 11 ^ ^ 15; this 

range corresponds to 1.6kpc <r < lOkpc and therefore contains 
red clump stars in the bulge and long bar. Their nature as an approx¬ 
imate standard candle, together with their abundance, provides the 
statistical power in the method. 

We model only the density and do not attempt to build dynam¬ 
ical models. In the future we will use the observational constraints 
from this work together with the made-to-measure method to con¬ 
struct self-consistent models whereby the stellar density generates 
the potential in which the stars orbit, similar to |Portail et al.|(20T5]). 
In principle fitting a non-parametric model similar to Wegg & Ger-| 
|hard| { |2013] ) would be preferable. Unfortunately the signal-to-noise 
of RCGs relative to the background of non-RCGs is significantly 
smaller in each field in the long bar than in the bulge region. The 
parametric models are fitted all fields simultaneously, thus improv¬ 
ing the bar signal. The planned made-to-measure modelling would 
again be non-parametric, but would connect different fields by re¬ 
quiring that the model be dynamically self-consistent. We regard 
the parametric models constructed in this section as a useful step 
in uncovering the properties of the data even though they do not 
statistically match the data perfectly. 

Our model density consists of three distinct parts: 


(i) An N-body model taken from |Portail et al.| ( |2015| l. These N- 
body models were adjusted to fit the density and kinematics in 
the galactic bulge at |/| < 10° and therefore fit this region well. 
As in [Portail et al. | ( |20 1 5] l the model is placed so that the bar is at 
an angle of a = 27°. 

We use an SPH cubic spline kernel to evaluate a smooth density 
field from the N-body model. For clarity we repeat the formulae 
here following the notation of [Hunt & Kawata| ( |2013| l. We use 
the cubic spline kernel 


W{s,h) 



1 — 6{s/6)'^ -\-6{s/Kf' 
2{\-s/hf 

0 


The density at a point Xi is then given by 


if 0 ^ s/h <1/2, 
if\/2^s/h< 1, 
otherwise. 

( 6 ) 


p{xi)= Y,miW{\xi-Xjlhj) 
i=i 


(7) 


where the sum runs over all N-body particles, mj is the mass of 
particle j. The smoothing length of particle j is evaluated from 
the local density through 


hj = ^ 



( 8 ) 


where T] is a parameter for which we have chosen rj = 3. Equa¬ 
tions]^ and are solved iteratively until the difference between 
iterations is less than 10“^. The choice of both these parameters 
mirrors [Hunt & Kawata| ( |2013| ) . 

(ii) We add parametric functions to represent the long bar compon¬ 
ent, the component of primary interest in this work. We choose 
the form (similar to |Freudenreich|19^|Picaud & Robin|2004[ 
with C|| = 1) 


Robin et al. 


2012 


P = 


^bar 

47rxoyozo 


exp 



- l/c±\ 

z 

jexp 

^0. 


X Cut 


^ ^out 


Cut 


Rin R 

f^in 


(9) 
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where x,y,z are right-handed galactocentric coordinates so 
that X is orientated along the bar major-axis, _y along the bar 
in termedi ate-axis, and z towards the north Galactic pole. R = 
is the galactocentric radius. We use a Gaussian cutoff 
for the inner and outer edges of the bar: 


Cut(x) 


exp(—ifx>l 
1 if X ^ 1 


( 10 ) 


Note that because of the Gaussian cutoffs the bar mass is less 
than the parameter M^ar- The inner cutoff is required because 
the N-body model already fits the magnitude distributions in the 
central region well as it was tailored to fit the data there using the 
made-to-measure method in |Portail et al.| ( |20T5] ). 

(iii) The Galactic disk scale length of the N-body model is 1.1 — 
1.2kpc, significantly shorter than that of the Milky Way, 

2 — 3 kpc (see |section A) . This is true in many N-body bar models 
formed from initially unstable disks where the bar length is typ¬ 
ically several disk scale lengths. That the long bar of the N-body 
model is less prominent than in the data is also related to the 
discrepancy in disk scale length, since moving outwards along 
the bar there are fewer stars than in the Milky Way. Here we 
do not attempt to model the disk to reconcile the discrepancy 
between the N-body disk and the Milky Way. Fitting the disk 
would be an involved task worthy of a separate work, while here 
we concentrate on features of the non-axisymmetric part of the 
density Le. the bar. Instead we add a best fitting additional expo¬ 
nential component to each modelled field magnitude distribution 
independently: 


^disk(MA:) =A'o,diskexp[^j(/iA:-13.5)] . (11) 

The result is that we are fitting the difference between the mag¬ 
nitude distribution and an exponential. We have verified that for 
the fields considered in this work, over the magnitude range con¬ 
sidered, that the Besangon Model ( [Robin et al.|2003| l predicts a 
nearly exponential magnitude distribution. Note that if the disk 
had a central hole that was significant for this work then the lack 
of RCGs in this region would result in a dip in the magnitude dis¬ 
tributions in a similar manner to the peak caused by the bar. Since 
this is not seen in the fields fitted here we consider [equation (fT)] 
adequate. 


In order to match the magnitude distributions in each field we 
convolve the resultant density from (i) and (ii) with a luminosity 
function to predict the number of stars as a function of jix before 
adding the disk from (iii) seperately. For clarity of notation we define 
the colour C = the corresponding reddening free colour 

Mc = {H- Ks) -E{H-Ks), the reddening Ec = E{H- TT,), and 
Rx = AkJE{H — Ks) which is a constant given by the extinction 
law. 

Consider <^>{MK,Mc)dMKdMc as the joint number of stars 
produced per unit mass between Mk to Mk + dMx and Me to 
Me + dMe where capital M denotes the absolute magnitude in the 
respective band. This two-dimensional colour-magnitude analogue 
to the luminosity function can be readily calculated by populating 
isochrones in a similar manner to the luminosity function. Define 
N{Ks,C) dKs dC as the number of stars observed in a pencil beam 
line of sight with solid angle ft) jointly in the range Kg to Kg + dKg 
and C to C + dC. Then the colour-magnitude version of the equation 
of stellar statistics is 

N{Ks,C) = Cojdrdp{r)'^[Ks -5\og{r/10pc) -AK,{r) , 

C-Ec{r)] (12) 


where Ak^ (r) and Ee{r) are the ^^-band extinction and the redden¬ 
ing inC = H — Kg along the line-of-sightrespectively. This can also 
be written in terms of distance modulus 

N{Ks,C) = j diils{ii)^[Ks-ii-AKXu) , C-Eciu)] (13) 

where A is given by A(/r) = (In 10/5) ft)pr^ expressed in terms of p. 
Changing variables from {Kg^C) to {jiK^Me) allows us to calculate 
N{iik) dfiK, the number of stars from to jj^k -{-dfiK- This is the 
quantity which we compare to the data. Using [elation (l)[ then 

N{ij^k) = JdMeN\^iiK-^RK (C-^c,rc) +^x,,rc , + 

= jduj dMcA{u}^[uK-U+RK{C-Mc^Rc) 

5 ^c] ■ ( 14 ) 

For an extinction law characterised by the constant Rk = ^{h-k ) 
the extinction can be calculated from the reddening. 


(r) = RKEe =Rk{C- Me) . (15) 


Therefore 

N{iik) = jdji J JMcA(p)d>[p/^-p+R^ (Mc-Mc rc) 

+ Af^,,RC , Me] (16) 

= j, (17) 

where (Af^u^) is the luminosity function in calculated from 

<^{MK,Me) through 

^idK= (dMe[M^ j^ -\-Rk {Me - Afc,Rc) + Af/^,,RC , Me] • 

(18) 

This differs slightly from the luminosity function in Mk because 
not all stars share the color of the red clump. Because we use the 
extinction free magnitudes (equation[^ extinction does not enter in 
the final equation provided the survey is complete and the extinction 
law is correct. This is true of all stars, not just RCGs. Practically 
the convolutions in [equation (17) [ are performed as parallel FFTs 
for speed. The exponential in N{fJ.K) used to represent the disk 
described in (iii) is then added to the result of [equation (17)[ 

The luminosity function is calculated from the BASTI 
isochrones ( [Pietrinfemi et al.|2004[ |. We fiducially use the lOGyr, 
a-enhanced isochrones together with a [Kroupa| p001[ l IMF and 
the metallicity distribution measured by [Zoccali et al. ( [2008[ l in 
Baade’s window to generate the luminosity function. To allow for 
the uncertainty in red clump magnitude, when fitted we allow a 
global shift in the luminosity function AK. 

We then fit by minimising the between the model and ob¬ 
served magnitude distributions assuming that the error is the Poisson 
error in observed number. If we label magnitude bins i and fields j 
then we minimise: 




2 


I I 

fieWs/ magnitu(iebins 7 


(V 


model /,7 ' 


NiJ 



(19) 


where A^odel/,; is the prediction of the model and Nij the observed 
number. We use bins in magnitude of 0.05 which are sufficiently 
narrow to not artificially broaden the luminosity function. We fit 
only fields where the exponential slope fitted to the background in 
[section 3| was greater than 0.55 since, as discussed in that section, 
smaller values indicate significant incompleteness at the faint end 
resulting in a smaller slope. 

Throughout we place the Sun a distance Rq = 8.3kpc from 
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the Galactic center ([Reid et al.[[201^ 

Chatzopoulos et al.[[2014[l 

and 25 pc above the Galactic plane (|A 

afz-Apellaniz|2001| Chenj 


|et al.|200T||Juric et al.|2008) , and Sgr A*in the physical Galactic 
mid-plane at Z? = —0.046°. This results in a slight tilt of the Z? = 0 
plane with respect to the physical Galactic mid-plane (Fig. 2 of 


[Goodman et al.|2014| ), but well within the errors when the Galactic 
coordinate system was originally defined. To move from axes v,y,z 
aligned with the principle axes of the bar to Galactic coordinates 
we first rotate by bar angle a about the z-axis, then move the sun to 
Rq = 8.3 kpc and zq = 25pc from the Galactic centre, and finally tilt 
the Galactic coordinate system by 0.12° towards Sgr A*( |Goodman| 
jet al.|2014l ). 

We describe the fitting as a three stage process: We first fit the 
N-body model to the central fields in jsubsection 5.1| we then add an 
additional component to fit the insufficient long bar component in the 
N-body model in jsubsection 5.2| and finally we add an additional 
component required to adequately fit the super-thin component 
towards the bar end in jsubsection 5.3| The model we regard as our 
best model is that with two parametric long bar components. We give 
the fitted parameters in [Table l] and the resultant physical quantities 
such as the bar length and the mass of the barred components in 
ITable 21 


5.1 Bulge Fitting 


To demonstrate and check the method we first fit the central bulge 
region. 


<10°. We take the model named M85 from Portail et al 


( j2015j l as a fiducial fit. We fit for two global parameters: the normal¬ 
isation of the model density G, and an offset in red clump magnitude 
from that assumed by the luminosity function constructed from 
the isochrones. 

We show in the top left panel of Fig. 12 the resultant value 


in each field across the entire range of longitudes, although only 
the central region was fitted. We also show in the top right panel 
of [Fig. 12] the mean absolute fractional error in each field over a 
magnitude range covering the red clump. Comparing these panels 
it can be seen although that the reduced x^ is formally poor in the 
central fields the fractional error is extremely small. The larger x^ 
here is due to the large numbers of stars and resultant small Poisson 
errors. The reverse is true at |Z?| > 8°: the small number of stars 
results in larger statistical errors and therefore larger mean absolute 
fractional deviations, however the x^ values demonstrate the fit is 
statistically good. It is evident that the model fits poorly in the fields 
close to the plane at / > 10°. This is a result of the insufficient long 
bar component in the N-body model in comparison to the Milky 
Way. 

In the left hand column of [Fig, l^ we show some example 
fields and the resultant fits. Again it is clear that the fit in the bulge 
region is excellent, however the long bar is insufficient outside this 
central region. The fact that the N-body model together with the 
disk model fits well in the bulge region shows that the disk model is 
reasonable, since the N-body model was constructed to fit bulge-only 
deconvolved data. 

We give our fitted parameters in [Table l] The best fitted shift 
in red clump magnitude, /S.K, is 0.17, suggesting that either Rq is 
slightly less than our adopted 8.3kpc, or bulge red clump stars are 
slightly brighter than assumed in the luminosity function. Since 
our assumed bulge luminosity function places the red clump at 
— 1.59, this would correspond to a fitted red clump absolute mag¬ 
nitude of M/^ rc = —1.75 — 5log[/?o/8.3kpc]. This is brighter than 
the Rc = —1.61 measured from nearby Hipparcos stars (Laney 
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Figure 13. Example fits to the magnitude distribution in selected fields. 
In the left column we fit only the N-body model (and the additional field- 
by-field disk component [subsection 5.1) . In the middle column we fit the 
N-body model together with a 1 component bar ( [subsection 5.2^ . In the 
right column we fit the N-body model together with a 2 component bar 
( [subsection 5.3) . The data with its associated Poisson error is plotted in 
black. The best fitting model is shown in blue. The model consists of a thin 
bar component (green), a super-thin bar component (orange), an exponential 
to represent the foreground disk (magenta) and the N-body model (grey). 


[et al.|201^ . The normalisation of the model density, G, of 1.47 
depends on the N-body model chosen which differ in their dark 
matter content in the Bulge, and is also degenerate with the IMF. 
This is discussed in detail in [Portail et al.| ( |20T5] ). 


5.2 One Component Long Bar Fitting 


We now add an additional parametrised function to represent the de¬ 
ficient long bar in the N-body model clear from jsubsection 5.1] The 
functional form given in jequation (9) [i s equivalent to the exponential 
shape function fitted in Robin et al. \20\2\. W e choose this form 
based on the results of isection 3l and lsection 4l that the vertical struc¬ 
ture is approximately exponential. The resultant horizontal density 
profiles are elliptical in the case of = 2 and boxy for > 2. We 
have added an additional Gaussian inner cutoff function since the 
N-body model already fits well in the central region. 

We fit this functional form, together with the N-body model 
and minimise x^ in fields with / > —10° . The best fitting parameters 
are given as the second row of |Table l] The reduced x^ nnd mean 
absolution deviation are plotted in the second row of [Fig. 12] from 
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Figure 12. In the left-hand three panels we show the reduced for the magnitude distribution in each field. In the right-hand three panels we show the mean 
absolute fraction error in each field over the range 13 ^ ^ 15. In the upper panel we show the fit using just the N-body model. In the middle panel we add an 

additional parametric bar component. In the lower panel add a second parametric component to fit the super-thin component at the bar end. White regions are 
those which were not fitted because they did not pass our completeness test. 


Table 1. Parameters found when fitting model densities to the magnitude distribution. The parameters and the method of fitting is described in | section's] Note 
that Mbar differs from the physical bar mass give in |Table 2| because of the Gaussian cutoffs. 





G 

logMbar/M© 

a (deg) xo(kpc) yo(kpc) 

zo(kpc) /?out(kpc) crout(kpc) /?in(kpc) 

CTin (kpc) 


N-body model only 

2.06 

0.17 

1.47 

- 

- 

- 

- 

- 

- 

- 

- 

N-body model + 1 component 

1.58 

0.17 

1.48 

10.15 

28.4 

3.05 

0.68 

0.09 3.85 

0.72 3.45 

0.71 

2.27 





f 10.87 

29.1 

0.82 

0.23 

0.20 4.37 

0.11 3.26 

0.56 

3.98 1 

N-body model -t- 2 components 

1.54 

0.16 

1.45 

\ 9.97 

30.0 

-32.66 

1.88 

0.04 3.36 

0.76 7.09 

2.24 

1.67 1 


Table 2. The derived physical parameters, bar mass 


and bar length, of the models with parameters fitted in |Table l] and described in |section3| 


Bar Mass (lO^M©) 


Bar Half Length (kpc) 



N-body 

Thin Bar 

Super-thin Bar 

Total 

'^drop 

L^m—l 

'^prof 

^mod 

N-body model only 

11.0 

- 

- 

11.0 

4.19 

4.03 

4.21 

- 

N-body model + 1 component 

11.1 

8.8 

- 

19.9 

5.50 

5.44 

5.65 

4.87 

N-body model -t- 2 components 

10.9 

4.0 

3.3 

18.1 

4.98 

4.96 

5.23 

4.73 


The bar mass of the N-body model is calculated by integrating the face-on surface density over all radii with the minor axis profile in the surface density 
subtracted. The thin and super-thin bar masses are the total mass in each component. The bar half lengths are as defined in |Athanassoula & MisiriotisH2QQ2t and 
all are measurement on the face on density map: Ldrop is the radius at which the ellipticity drops fastest, is the radius at which the m = 2 component of the 
face on image drops below 20% of its maximum, and Lpj-of is the radius beyond which the major and minor axis density profiles agree within 30%. Our threshold 
for Lprof is larger than the 5% used by|Athanassoula & Misiriotis|J2002| which we found gave spuriously large bar lengths for the densities considered here, 
f-mod is the point at which the difference between the major and minor axis face-on surface densities falls to below 1/e of a fitted exponential long bar surface 
density profile. It has the advantage that it is independent of the axisymmetric disk. This measurement is spurious for the pure N-body model since it does not 
occur within the solar radius. 
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which it is clear that the resultant fit is greatly improved close to 
the plane in the region 10° < / < 30°. Note that the region near 
5° < / < 15° and \b\ ^ 5° is not yet well fitted. The fitted bar angle 
is a = 28.4° in agreement with the bar angle in the bulge region 
found by Wegg & Gerhard (20131 of a = (27 ±2)°. 

The additional mass associated with this component is 8.8 x 
\0^Mq for our luminosity function calculated from the BASTI 
isochrones for a lOGyr population with a |Kroupa|p001| l IMF . If 
we also include the non-axisymmetric mass from the N-body model, 
calculated by integrating the face-on surface density over all radii 
with the face-on minor axis profile subtracted, then the total bar 
mass is 1.99 x lO^^M©. Note though that the mass derived here, 
which is constrained by the number of red clump stars, is degenerate 
with the IMF and population age. This is because the mass per 
red elump star varies in a similar manner to the mass-to-light ratio. 
Refitting changing from our fiducial Kroupa IMF to a Salpeter 
IMF ( |Salpeter|1955| l increase all masses by a factor of 1.43, while 
keeping all other parameters the same. Changing from a lOGyr to 
a constant star formation rate would reduce the mass by a factor 
of 2.0. In this case since the luminosity function changes slightly 
the best fitting parameters change marginally from our fiducial 
values, however the physical properties of bar angle and length do 
not change significantly.Uncertainty in IMF and population age 
dominate the bar mass measured in this work and are larger than the 
differences between models or bar parameterisations. 


5.3 Two Component Fitting 

The adopted parametric density in |equation (9)| has a single constant 
scale height, however the results of |section^| suggest that the seale 
height is not constant becoming thinner towards the bar end. For 
this reason we have added a second component parameterised with 
the same functional form to represent the super-thin component. 
We construct the luminosity function of the super-thin component 
assuming a constant star formation rate, as opposed to the lOGyr 
age of the bulge and 180 pc scale height thin bar. As discussed in 
the previous section this choice reduces the mass-per-clump star 
compared to a simple lOGyr population. 

We initially allowed all parameters of the fit to vary, however 
we found that this resulted in an unsatisfactory fit. The higher num¬ 
ber of stars in the fields nearer to the bulge region result in smaller 
Poisson errors and therefore higher weight placed on fitting these 
well, so that outside this region the fit is less good. In particular the 
thin bar component became too thick in order to fit the fields near 
the bulge, and the resultant vertical structure in plots similar to Figs. 
[8}fTT]was poorly represented. 

Instead, motivated by |Fig. 9| we fixed the vertical scale height 
of the two components, finding that exponential scale heights of 
200 pc and 40 pc produced a reasonable representation of the vertical 
structure. The resultant fit is considerably improved, particularly 
near the bar end and between 10° < / < 20°, as can be seen in the 
lower panels of [Fig. 12[ We also observe this model similarly to 
the data, in slices from above, and plot the result as the lower half 
of |Fig. 5 1 finding it does a remarkably good job of reproducing the 
features in the data. We show in |Fig. 14] the fitted two-component 
model, separated by component, projected both from above and 
perpendieular to the bar along the minor axis. 

The angles of the component are a = 29.1° for the thin com¬ 
ponent and a = 30.0° for the super-thin component, again consistent 
with measurements of the bar angle in the bulge region. The mass 
of the 200pc thin component is 3.3 x IO^Mq while the mass of 


the 40pe super-thin component is 4.0 x lO^M© giving a total mass 
together with the N-body bar of 1.81 x lO^^M©. 


5.4 Bar Angle 


In our fiducial model the fitted bar angle of the parametric long bar 
ranges from a = (29 — 30)° depending on the number of compon¬ 
ents ( Table l| l. Fitting the other N-body models described in |Portail| 
|et al. “ pOlSI l with the same procedure gives angles in the range 


a = (28 — 31)°. Altering the age and metallically of the bar com¬ 
ponents while holding the bulge luminosity funetion fixed ehanges 
the fitted bar angle slightly. This is because the RCG luminosity 
then changes between the inner galaxy and the long bar region. Mak¬ 
ing the super-thin component super-solar metallicity ([z/H] =0.17) 
increases the fitted angle of this component to 32.4°, while chan¬ 
ging the 200pe component to have a constant star formation history 
increases its angle to 32.7°. We therefore increase our range of 
possible long bar angles to a = (28 — 33)° to eneompass this. This 
is consistent with alignment with the angle found in the central 
|/| < 10° in Wegg & Gerhard| ( [2013| l of a = (27 zb 2)° which is our 
assumed N-body bulge angle, and with a = (29 — 32)° found by 
|Cao et al.1p013| ) for a simpler parametric |Dwek et al.| ( |1995] ) bulge 
model. 

We regard the one component long bar model described in 
[subsection 5.2| as a useful confirmation that the fitted bar angle is 
insensitive to the parametric model. This is a significantly simpler 
model that fits the data less well than the two component long bar 
described in [sub section 5.3 [ but still recovers a consistent bar angle. 

We have performed an MCMC to estimate errors on both the 
model parameters, and the physical properties such as bar mass 
and length. The resultant statistical errors are extremely small, sig¬ 
nificantly smaller than the difference between fits with different 
models. We therefore consider the statistieal errors to be negligible 
in comparison to systematic errors quoted above. We have also 
refitted starting from different initial conditions to verify that the 
fitted parameters such as the bar angle is not a local minimum, or 
strongly dependent on initialisation position. While the same min¬ 
imum in was not always reached by the minimisation procedure 
the difference was smaller than between models. 


5.5 Bar Length 

In external galaxies and N-body simulations the bar length is depend¬ 
ent on definition with many possible (g.g. jAthanassoula & Misiriotis] 
|2002| ). Even in N-body simulations with strong bars, where the full 
6-dimensional phase space of every particle as a function of time 
is available with negligible error, the different definitions produce 
different bar lengths, typically at the level of ^ 10% for strong 
bars. We therefore do not expect an unambiguous and definitive bar 
length for the Milky Way, particularly with the more challenging 
viewing geometry. Instead we estimate the bar length of our best 
fitting model densities using several definitions. These bar length 
measurements are summarised in ITable 2l 

Specifically we use three bar length definitions found in jAth-j 
janassoula & Misiriotisj ( |2002j l to give reasonable bar lengths of 
N-body models: (i) the point where the elliptieity of the face 
on profile drops fastest, (ii) L^= 2 , the point where the m = 2 Fourier 
mode of the face on image drops to 20% of its maximum, (iii) Lpj-of 
the point where the major and minor axis profiles agree within 30% 
(this a larger threshold than the 5% used by jAthanassoula & Misiri-j 
|otis|2002| whieh we found gave spurious results for our densities). 
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Fitted Model N-body model + Besangon Disk + Thin Bar Component + Super-Thin Component 



Figure 14. The three upper panels show the view from the north Galactic pole of the surface mass density of the best fitting two component model described in 
|5.3| The lower panels show the same model observed side-on. We also add the Besan^on model disk density described in |Robin et al.| ( |200^ since the disk of 
the N-body model alone is insufficient. The left hand panels show the N-body model together with the Besangon model disk. The central panels additionally 
include the 200 pc thin bar component. The right hand panels additionally include the 40 pc super-thin bar component. The bar length was measured from these 
images using the methods described in |subsection 5.5| 


We add one additional bar length measurement, Ljnod- We take 
the difference between the face-on major and minor axis surface 
density profiles along the bar and fit an exponential to the long 
bar region. We then define the bar length as the point at which the 
density falls to \/e of the exponential profile. In the case of an 
analytic exponential bar with a Gaussian cutoff, like our parametric 
long bar functions, this corresponds to defining the bar length as 
Ro~\~ 


All these methods were applied to face-on images of our densit¬ 
ies to which we added the density of the disk in the Besangon model. 
We show these face-on images in |Fig. 14] and the resultant bar half 
lengths are given in [Table '2\ All our stated bar lengths are the half 
length, defined as the distance from the galactic centre to the bar 
end. 


Before considering the bar length measurement we first return 
to the data near the bar end. In |Fig. 15] we show histograms stacked 
in Galactic latitude as a function of longitude. We show both the 
positive and negative longitude sides to demonstrate the peak at 
positive longitudes is non-axisymmetric. To make the plot clearer 
we also subtract an exponential in f^K which can be thought of as 
representing the background of non-RCGs. The bar is clear and well 
localised to / < 26°. At / > 30° while non-axisymmetry still appears 
it is much less significant and fainter than would be expected for 
the bar. In the region 26° < / < 30° the non-axisymmetric excess 
weakens, broadens and becomes fainter. We therefore presume that 
the bar ends in this region, possibly transitioning into the spiral 
arms. If we convert these longitudes of the bar end to a bar length 
assuming that the bar lies at a ^ 27° and that projection effects 
are negligible we would recover a bar half length between 4.4 and 
4.8 kpc. 

We produced a similar plot to |Fig. 15| showing just the N-body 
model and it is clear that the bar ceases to be significant at too low 
longitude. In contrast plotting the one component long bar model it 
is clear that the bar extends beyond the data in longitude to where 


there is no non-axisymmetry in the data. For this reason we disregard 
the bar length measurements of these models in [Table 2] 

Instead the two component model is a significantly better fit to 
the stacked data near the bar end in [Fig. 15] We show it compared 
to the data together with variations in which the bar was artificially 
lengthened and shorted by adjusting the outer cutoff by 0.5 kpc. The 
model with an artificially shortened bar is insufficient particularly 
beyond / = 25°. In contrast the model with an artificially lengthened 
bar predicts excessive non-axisymmetry beyond / = 30° when the 
positive and negative latitudes have similar counts at the distance of 
the bar. 

Because the two component model appears to reasonably fit 
the stacked data in [Fig. 15] we consider the measurements of this 
model to be our fiducial bar half length. These measurements lie 
in the range 4.73 — 5.23kpc. We have repeated this process on the 
other N-body models finding that the variation between models 
is smaller than the variation between methods of measuring bar 
length. The one component bar length appears longer, however it is 
evident from [Fig. T^ that this model fits poorly in the region beyond 
/ > 30°. Therefore taking the average and standard deviation of 
these measurements we consider our fiducial estimated bar half 
length for the Milky Way to be (5.0 ± 0.2) kpc. 


6 DISCUSSION 


6.1 Continuity Between Box/Peanut Bulge and Long Bar 

Two lines of evidence in this work support that bar and bulge appear 
to be naturally connected: the angle between the Box/Peanut (B/P) 
Bulge and Long Bar is small, and the scale height along the bar 
decreases smoothly. 

We find in this work that the long bar has bar angle in the 
range a = (28 — 33)° consistent with recent determinations of angle 
found at |/| < 10° in the B/P Bulge (e.g. Wegg & Gerhard|[2013 1. 
We find this angle by fitting the magnitude distributions through 
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Figure 15. Histograms of stars near the bar end. In black we show histograms of stars with positive longitudes and 0.15° ^ \b\ < 1.35° in one degree longitude 
intervals. The histograms of jik jequation (1)} were converted to distance assuming all stars are red clump stars and an exponential in jik was subtracted for 
visualisation, which can be thought of as representing non-red clump stars. In orange we show the same procedure applied to the data at negative longitues to 
highlight the non-axisymmetric features. In red we show our best fitting density model from |section 5| while in blue and green we show the same model but with 
the bar length reduced and increased by 0.5 kpc respectively. 


the parametric density functions in section 5| It was already clear 
however directly from the data in |Fig. 5 [ that the bar angle between 
10° < / < 20° was close to the a ^ 27° found in the B/P Bulge. 

The alignment of the B/P Bulge and long bar in this work 
is in contrast to some previous claims of misalignme nt (|Lopez 


Corredoira et al.|2006l|Cabrera-Lavers et al.|2007||2008| t. [Benjamin 


et al. 


( [2005 | l found a large angle of the long bar of a = (44 ±10)^ 


using GLIMPSE data. However the GLIMPSE long bar data was 
subsequently analysed finding a = (38±6)° (Zasowski|2012| or 
a = 35°, jZasows^ et al.|[20T^ due to a fainter assumed RCG 
absolute magnitude. This results in the derived distances to the 
long bar being reduced and therefore the bar angle decreasing for 
the same Sun-Galactic center distance, Rq. Note that because we 
use data that extends in longitude to the bulge region and allow 
the red clump absolute magnitude to vary in our fitting we are 
insensitive to this degeneracy in this work. Some misalignment 
claims appear to be a result of using the endpoints to derive the 
bar angle {e.g. [Gonzalez-Eemandez et al.|2Q12[|de Amores et al.j 
|2013j ), however the projection effects can be extreme especially for 
the far end of the bar. Large bar angles from fitting the distance 
of clump giants (g.g. jCabrera-Lavers et al.|2007||2008] l are more 
difficult to reconcile with this work and may be related to the very 
different selection criteria. In those works RCGs are selected from 
the colour-magnitude diagram together with an extinction model 
which predicts the reddening as a function of distance modulus, 
while in this work they are statistically identified as an excess above 
the smooth background of non-RCGs. 

In addition we find that the scale height shown in jFig. 9| appears 
to smoothly transition between the B/P bulge and the long bar which 
also suggests that both are part of the same connected structure. This 
transition is similar to N-body models. To demonstrate this we show 
in jFig. 16| the scale height of one of the initial barred N-body models. 
The transition in vertical structure from a short central scale height, 
to a large scale height through the 3D B/P region, to a short scale 
height in the long bar can be seen to arise naturally in the N-body 
model prior to fitting the bulge data. 

Given the near alignment it seems natural that the long bar is the 



Figure 16. The scale height of the N-body model M85 from |Portail et al.| 
( |2Q15) as a function of distance along the major axis of the bar. We use a 
secir scale height since this is a better fit to the N-body vertical structure 
than an exponential. The model shown is prior to fitting the bulge data to 
demonstrate that the transition from short central scale height, to large scale 
height in the B/P region, to short scale height in the long bar region arises 
naturally. 


in-plane extension of a vertically extended inner part of the bar, the 
Box/Peanut bulge. This is similar to the structures found in N-body 
simulations where the buckled three-dimensional boxy/peanut bulge 
is shorter than the bar {e.g. |Athanassoula|2005||Martinez-Valpuesta 
et al. 2006| l, as well as in external galaxies ( Erwin & Debattista 
20131. This situation was previously suggested in the Milky Way 
by [Martinez- Valpuesta & Gerhard] ( |20TT] l and [Romero-G 6 mez et al.j 
( [201 1| despite the possible misalignment. 

In the lower panel of jEig. 14] we show the side-on projection 
of the best fitting model. These side projections where the three 
dimensional peanut has a lesser extent than the bar is similar to 
the side projections in [Athanassoula[ ( [2005] l and the observations of 
[Bureau et al.| ( [2006] ). 
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6.2 Thin and Super-thin Bar Component 


From fits to the vertical profile of the RCGs we find evidence for 
two components associated with the long bar of the Milky Way. A 
component with a scale height of ^ 180pc which we term the ‘thin 
bar’ since it would appear to be formed from the counterpart of old 
solar neighbourhood stars, and a second ‘super-thin bar’ component 
with scale height 45 pc. 

The density better fitting a broken exponential alone would 
be weak evidence for the presence of two distinct components: As 
demonstrated by the controversial thin/thick disk discussion in the 
solar neighbourhood it is possible for a continuum stellar distribution 
to mimic a double exponential in vertical profile {e.g. |Bovy et ^ 
|2012| t. However the relative contributions of the thin and super-thin 
components change along the bar. The thin component decreases 
outwards approximately exponentially with a scale length similar 
to the Milky Way disk, while the super-thin component increases 
outwards towards the bar end. This spacial mismatch strengthens 
the argument for two distinct components. Additional abundance 
and age measurements would strengthen the argument further. 

The short scale height for the super-thin bar is similar to the 60- 
80pc found in X^-band imaging of NGC 891 by |Schechtman-Roo^ 
|& Bersha'd^ ( |2013| l. It is slightly lower than the expected value in 
the solar neighbourhood. Young stars locally have a scale height 
of ^ 60pc {e.g. |Joshi|2007| from measurements of OB stars and 
young open clusters) and the vertical disk heating in the last 2Gyr 
is insignificant |Holmberg et al.| ( |2007| ). 


If we extrapolate the local thin disk properties inwards we 
can predict the resultant velocity dispersion of the components. 
Assuming that the thin bar’s scale height is set by its self gravity 
then the vertical dispersion, cj^, will vary as oc Z/i where Z is the 
surface density and h the scale height. Assuming that the surface 
density Z is exponential Z oc exp {R/Rd) then the surface density at 
Galactocentric radius 4kpc is a factor 7.4 higher than in the solar 
neighbourhood for Rd = 2.15kpc l |Bovy & Rix|2013|l. Si nce the 
local thin disk scale height {e.g. 225pc: |Veltz et al.|2008| 300pc 
|Juric et al.|2008| l is slightly larger than the thin bar scale height the 
resultant dispersion would then be a factor ~ s/JA x 180/225 = 2.2 
larger, or approximately ~ 4 6kms~^ assuming a local vertical 
velocity dispersion of 21 kms“^ { [Binney et al. 2014| ). 

Estimating the velocity dispersion of the super-thin component 
is less certain since it depends on the extent to which its vertical 
structure is governed by self-gravity. We consider the two extremes 
of completely self-gravitating and non-self-gravitating and expect 
these to bracket the true structure. Assuming that it is self-gravitating 
a similar extrapolation gives ~ 30km s~ ^, assuming locally that 
stars younger than 2Gyr have ~ 10km s“^ jHolmberg et al. 2007| t 
and a scale height ~ 60pc ( |Joshi|2007] t. The dispersion of a non-self- 
gravitating structure would be lower. Assuming that the potential 
is near harmonic 0 oc vertical frequency is set by 

the thin component then the vertical frequency will be oc cj^//? 
and hence 2.7 times higher than in the solar neighbourhood. The 
resultant dispersion for the super-thin component would be (7^ oc y^h 
and therefore ^ 2.7 x 50pc/60pc x 10kms“^ 22kms“^. 

The origin of the super-thin component is unclear. The small 
scale height means the stars cannot have experienced much vertical 
heating and therefore should be younger than the thin component. 
To have formed RCGs the stars must still be > 500 Myr old but 
galaxies with ongoing star formation have a strong bias to younger 
RCGs with a peak of age ^ 1 Gyr ( [Girardi & Salaris|2001[|Salan^ 
|& Girardi|2002|l. They cou ld be related to star formation towards 
the bar end ( [Phillips|1996| l. Because their lifetime is longer than 



Figure 17. Our fitted two component Milky Way model density ( |subsec-| 
|tion 5.3) observed as if an external galaxy with bar angle 45° and inclined at 
i — 60°. The line-of-nodes of the disk is horizontal. As described in lErwinI 
|& Debattista| ( |2Q13') , this is an orientation where both the boxy and barred 
zones are visible with the angle of the boxy zone closer to the line-of-nodes 
than the projected angle of the long bar. The full length of the bar and boxy 
zone were measured from this image visually in a manner similar to |Erwin| 
|& DebattistaH2013} and are marked in blue and red respectively. 


the orbital timescale the interpretation would be that they are stars 
on bar following orbits that formed towards the bar end and spend 
more time towards their apocenter near the bar end. Even without 
any star formation, since the bar grows with time (e.g. [Martinez^ 
IValpuesta et al.|2006] l, the bar will have grown into the star forming 
disk. These recently captured bar stars would then spend more of 
their orbital period at their apocenter near the bar end and could 
therefore have a similar density profile to the super-thin component. 
Eurther dynamical and chemical information is needed to distinguish 
between these possible scenarios. 


6.3 Bar Length 

Our estimated bar half length of (5.0 ± 0.2) kpc is larger than most 
previous estimates (e.g. 4.5kpc by |Cabrera-Lavers et al.|2008) . This 
is partly because our bar angle is less than most previous estimates. If 
one assumes that the super-thin component is not part of the bar the 
same bar length estimators applied only to the thin component would 
lead to an only slightly shorter bar length of (4.6 ± 0.3) kpc. Without 
dynamical information the exact nature of this region near the bar 
end and the possible transition to spiral arms or ring structures is 
difficult to determine. 

A further uncertainty is that in estimating bar lengths we use the 
face on view of our model together with the Besangon model disk. 
If the Besangon disk density is not accurate in this region it would 
slightly impact our estimated bar length. Eor example increasing the 
density by a factor of 2 reduces the bar length from (5.0 ± 0.2) kpc 
to (4.6 ± 0.2) kpc. This is why Ihe point where the difference 
between the major and minor axis profiles of the face on surface 
density drops to Xje oi the otherwise smooth and slower decline 
along the bar, is useful. Since it considers only the non-axisymmetric 
part of the density it is independent of details of the disk. 

A firm lower limit to the bar length can be determined visually 
from |Eig. 3| it is clear from the slices at Z? > 2° that the bar extends 
to / > 22° and therefore given a bar angle a = 30° the bar cannot 
be shorter than 4.0kpc. 

External galaxies show a strong correlation between the size of 
their boxy bulge, and their in-plane long bar which is seen as ‘spurs’ 
extending from the boxy region. [Erwin & Debattista| ( |2013| l finds 
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the ratio of the size of the boxy 3D bulge to the thinner in-plane 
bar is 0.42 zb 0.07. This was estimated by measuring a sample of 
moderately inclined external galaxies. We therefore inclined our 
model density and measured the boxy zone and spurs in a similar 
manner. In |Fig. 17] we place an external observer of the Milky Way 
inclined at 60° with the bar at APA = 45° to their line of sight. We 
find they would measure a box bulge size of 2/?box = 2.3 kpc, a bar 
length 2 Lbar = 8 . 8 kpc and therefore a ratio to the bar length of 0.26. 
This is towards the lower end of the external galaxies measured by 
[Erwin & Debattist^ ( |2013| ), at the edge of the observed range. 

The bar length measured in this work here has implications for 
the pattern speed of the bar. The dimensionless bar rotation velocity 
is ^ where Rqr is the corotation radius. Since the bar can¬ 

not extend beyond corotation ^ ^ 1 and our measurement of bar 
length implies corotation must lie outside (5.0 zb 0.2) kpc. In turn this 
limits the pattern speed to be lower than < (45 zb 2) kms“^ kpc“^ 
for a flat rotatio n curve near the b ar end with circular velocity 
Vc = 218kms“^ ( Bovy et al.|2012 >. This is in some tension with 
the pattern speed derived form the interpretation of the Hercules 
stream as being due to the outer Linblad resonance of the bar. jAnt-j 


oja et al. ( [2014| l find a pattern speed of (48.2 zb 0.5) kms ^ if the 


bar lies at a = 29° with Rq = 8.3kpc and circular velocity in the 
solar neighbourhood of = 220kms“^. Reconciling these results 
requires either an extremely fast bar with ^ ^ 1 , where the bar 
ends at corotation, or that the Hercules stream has a different ori¬ 
gin than the outer Lindblad resonance. In contrast, by constructing 
dynamical models of the B/P bulge, jPortail et al.j ( |2015| l found 
25 — 30km s“^ kpc“^. Combining this with our bar length and its 
error gives ^ in the range 1.4 — 1.9 and therefore the Milky Way 
would be a slow rotator (^ > 1.4, |Rautiainen et al.|2008| >. 

In the future we plan to combine the constraints from this 
work with kinematic data in the long bar region and perform made- 
to-measure modelling in a manner similar to jPortail et al.| ( |20T5] l. 
Constructing self-consistent models of the long bar with added 
kinematic data will better indicate which stars are on bar following 
orbits, and thus constrain the dynamical structure near the bar end 
and the bar length. 

In addition the Gaia mission i jde Bruijne|20I2| | is capable of 
some elucidation of the long bar. Extinction will limit the effective¬ 
ness of its broadband visual measurements near the galactic plane. 
However the result that the long bar extends out of the plane means 
that at |Z?| > 1 ° it should be possible to measure reliable parallaxes 
and proper motions for a sufficient number of bright RGB stars. In 
addition GAIA will allow for better characterisation of RCGs as 
standard candles, and therefore further their use in studying Galactic 
structure in the higher extinction regions, where NIR measurements 
are required. In particular, better characterisation of the population 
effects on RCGs will be of importance in the future since, for ex¬ 
ample, our largest uncertainty in bar angle is due to the uncertainty 
in the difference in RCG magnitude between bar components. 


models that simultaneously fit all fields we reach the following key 
conclusions: 

• The bar extends to / 25° at |Z?| 5° from the Galactic plane, and 

to / ~ 30° at lower latitudes. 

• The long bar of the Milky Way lies at an angle of a = (28 — 33)°. 
This is consistent with being aligned to the Milky Way’s Bulge (e.g. 
I Wegg & Gerhard|2013[|Cao et al.|2013] l. We find this angle from 
the parametric density fits, but it is also visually evident from the 
raw data in jFig. 5| 

• The overall scale height of the bulge transitions from a short central 
scale height, to large scale height in the B/P region, to short scale 
height in the long bar region. This transition arrises naturally in 
N-body models of barred galaxies, and together with the alignment 
of the long bar and B/P region indicates that the bar and bulge are 
part of a single innately connected structure. 

• We find evidence for two scale height components. We find a 
^ 180pc ‘thin bar’ component, which decreases in density out¬ 
wards, and which we interpret as the barred counterpart of the solar 
neighbourhood thin disk. In addition there is a ^ 45pc ‘super-thin 
bar’ component whose density increases outwards along the bar. 
This component is likely to be related to younger stars (~ 1 Gyr) 
towards the bar ends. 

• The offset in b of the bar is consistent with symmetry about the 
physical Galactic mid-plane under the assumption that the sun 
lies 25 pc above the Galactic plane and Sgr A*lies in the phys¬ 
ical Galactic plane. The vertical structure is also consistent with the 
bar lying in the physical Galactic mid-plane to < 5 pc. 

• We find a bar half length of (5.0 zb 0.2) kpc measured by applying 

commonly used bar length definitions to our best fitting parametric 
model. This agrees with simple estimates directly from the data near 
the bar end ( [Fig. 15| l. The bar length is still somewhat ambiguous 
since it depends on the interpretation of the nature of the super-thin 
bar, and excluding this results in a bar length of (4.6 zb 0.3) kpc. It 
is also visually clear from the data in jFig. 5 [ that the bar appears 
straight extending to at least / 22° and therefore cannot be shorter 

than 4.0kpc as a lower limit. 

• Our estimated total bar mass in our fiducial parametric model is 
1.8 X 10^^Mq. This arises from the number of clump stars in the bar 
together with a ‘mass-to-clump’ ratio arising from our assumed IMF 
and component ages. Varying these assumptions dominates the error. 
For example changing from our assumed jKroupa j ( |200 1 j l IMF to a 
jSalpeterj l jI955j l IMF increases the long bar mass by a factor 1.43, 
while changing from an old 10 Gyr stellar population to a constant 
star formation rate reduces the mass by a factor of 2.0. 
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7 CONCLUSION 

We have investigated the structure of the Milky Way’s bar outside 
the bulge, termed the long bar, using red clump stars as a standard 
candle and tracer of the underlying density. We have combined data 
from UKIDSS, VVV, 2MASS and GLIMPSE and have constructed 
magnitude distributions in many fields covering the central |/| < 
40°, \b\ <9°. We concentrated on the ^^-band and corrected for 
extinction star-by-star. On the basis both of fitting the clump stars 
in each field individually, and of constructing parametric density 
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APPENDIX A: TRANSFORMATIONS TO 2MASS 
PHOTOMETRIC SYSTEM 

Al Transformation from UKIDSS 

The UKIDSS survey is performed in the WFCAM photometric 
system. From DR2 onwards the zero points were calibrated from 
2MASS but the photometric system is slightly different (equations 
4-8 of [Hodgkin et al.|2009| l. The most straightforward method of 
transforming WFCAM magnitudes to 2MASS would be to use the 
inverse of the calibration transformations. However these rely on the 
/-band magnitude and the /-band would then be the limiting band 
in high extinction regions, having almost twice as much extinction 
as the //-band. Instead for the H and K bands we measure the 
transformation from the // — /f color for giants, which can then be 
applied directly in equation 

To derive the transformation we cross match 2MASS and 
UKIDSS sources within F'and Imag with 2MASS photometric 
quality classified as ‘A’ or ‘B’ in the /, H and K bands. A cut was 
applied in color-color space of UKIDSS to remove nearby dwarfs 
and select mostly giants: 0.25 < {J — K) — 2.5(// — K) < 0.55. In 
order to avoid saturation in UKIDSS and Eddington bias (over rep¬ 
resentation of the more numerous faint objects due to measurement 
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Figure Al. Comparison of photometry from 2MASS (subscript 2) and 
UKIDSS (subscript W) for the K and //-bands. Contours show the density 
of cross matched sources selected as described in the text. Red points are the 
binned sigma clipped mean and the orange line is the derived photometric 
transformation. 



Figure A2. Comparison of photometry from 2MASS (subscript 2) and 
UKIDSS (subscript W) for the /-band. Contours show the density of cross 
matched sources selected as described in the text. Red points are the binned 
sigma clipped mean and the orange line is the derived photometric trans¬ 
formation. 


error) in 2MASS we also select only sources with l\ <K < 13.5, 
12.5 < // < 15 and 12 < / < 15.5. 

With this sample of cross-matched sources we then derive the 
photometric transformations. We divide the sample into 20 color 
bins each with equal numbers of stars and in each calculate the 
difference in magnitude between 2MASS and UKIDSS using an 
iterative sigma-clipped mean at the 2.5 sigma level. Our derived 
transformation is then the linear regression of these points, neglect¬ 
ing the first and last bins to reduce outliers. This process is shown 
applied to the H and K bands in figure |Al| and the / band in |A2| 


The resultant transformations are: 

Ji — Jw ~ 0.02(/;y — //w) ~ 0.03 
H2=Hw + ()-\2{Hw-Kw) (Al) 

K2 = Kw + 0.02(//vy - /fvp) + 0.01 

where subscript 2 refers to the 2MASS system and W the WFCAM 
system used in UKIDSS. 

A2 VVV Transformation 

The VVV survey is performed in the VISTAWIRCAM photometric 
system. As for the WFCAM system this is tied to, but different 
from the 2MASS systenQ In principle the same method as for the 
WFCAM system could also be used for VISTA. However the VVV 
system is significantly closer to the 2MASS system making this less 
necessary. There is however a field-to-field scatter in high extinction 
regions where there can be variations in zero point of ~ 0.1 mag. We 
instead take a different approach and as in |Gonzalez et al.|p01 l| t and 
|Wegg & Gerhard] ( [201 3 | t we re-estimate the zero points for each field 
by cross matching bright but unsaturated VVV stars with 2MASS. 


APPENDIX B: GLIMPSE DATA 

In this section we show the results of GLIMPSE 3.6qm and 
4.5 qm data with a similar analysis as applied to the K-band in 
the main text. Although the GLIMPSE data does not have the wide 
coverage of the K^-band data, the reduced extinction in GLIMPSE 
makes it useful to corroborate the results of the NIR data. In |Fig. B l] 
we plot the surface density of stars over a narrow range of extinction 
corrected magnitudes, the equivalent to |Fig. 1| In |Fig. B2| we show 
the fitted parameters of the statistically identified red clump stars: 
their number density Arc, the distance Jrc calculated from and 
their dispersion oi-c- This figure is the equivalent to |Fig. 4| in the Kg- 
band. Completeness limits the reliable identification of RCG stars to 
/ > 8°. In |Fig. B3] we show the equivalent top down view of the raw 
data as was made for the K^-band in |Fig. 5| As for the K^-band a bar 
angle of a = 27° appears close to the data while 45° does not. Near 
the dispersion of the fitted RCG is larger than the K^-band. This is 
expected because the RCG fits lie close to the magnitude limit of 
GLIMPSE in this more crowded and higher extinction area. 


^ For details see http://casu.ast .cam.ac.uk/surveys-projects/ 
vista/technical/photometric-properties 
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Figure Bl. The surface density of stars in the GLIMPSE [3.6jU]-band over the extinction corrected magnitude range 10.5 < [3.6/i]o <11. Asymmetric number 
counts in I close to the plane demonstrate non-axisymmetry. Extinction is corrected using the RJCE method i.e. the H — [4.5jU] colour excess as in equation]^ 
and data outside the colour bar range are plotted at its limit. Grey areas are those not covered by the GLIMPSE survey. The equivalent plot for the ^-band is 
shown in Fig.[^ 
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Figure B2. As Eig.j^but for the GLIMPSE 3.6p.m (left hand panels) and 4.5p.m data (right hand panels). 
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Figure B3. As jEig. 5| for the long bar of the Milky Way viewed from above using the GLIMPSE data in the 3.6p.m band (above) and 4.5p.m band (below). 











































































































